arXiv:1503.01504vl [math.ST] 5 Mar 2015 


On uniform performances of random 
polytopes and their functionals in convex 

bodies 


Victor-Emmanuel Brunei 

Cowles Foundation for Economics, Yale University 

victor-emmanuel.brunel@yale.edu 


Abstract 

In a previous paper we considered the random polytope defined as 
the convex hull of n i.i.d. random points uniformly distributed in a 
convex body and we focused on the Nikodym distance between that 
convex body and the random polytope. We proved a uniform deviation 
inequality, yielding tight moment inequalities in a minimax sense. In 
this paper, we use a different technique in order to prove similar types 
of inequalities for the Hausdorff distance between the random polytope 
and the convex body. We no longer restrict ourselves to the case of 
uniform distributions. However, our results allow to recover and im¬ 
prove some known inequalities for the convex hull of uniform points in 
a convex body or on its boundary. We also prove moment inequalities 
for some class of functionals of the random polytope, including mean 
width as a special case. 

Keywords, convex body, convex hull, deviation inequality, Hausdorff dis¬ 
tance, random polytope 

1 Preliminaries and notation 

1.1 Introduction 

Among all types of random polytopes defined in the litterature, we focus 
in this paper on convex hulls of finitely many independent and identically 
(i.i.d.) random points whose probability measure is supported on a convex 
body K in the Euclidean space M d . From now on, we denote this convex 
hull by K n , where n is the number of random points. We often refer to 
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these points as observations and to K as their support. The dimension d 
of the ambiant space is a fixed positive integer. We call convex body in 
M. d any convex and compact subset of M. d with positive volume, i.e., of full 
dimension. We are interested in the performance of K n as an estimator 
of the body K. In [6], we considered the case where uniformly distributed 
observations in K and we studied the Nikodym distance between K n and 
K, i.e., the volume of the symmetric difference between K n and K. We 
proved the existence of two positive constants r and c which depend on the 
dimension d only, for which the following holds. For all convex body K c W d , 
all n > 1 and all positive real number x, 


K I x 

— < rn d+i h — i 
n 

_c_ 

with probability at least 1 - ce d d . Here, | • | stands for the volume and A is 
the symmetric difference. Note that the denominator d d in the exponential 
term becomes very large when the dimension d grows and we believe it could 
be much improved by refining the technique of our proof. 

In the present work, we investigate the behaviour of the Hausdorff dis¬ 
tance between K n and K, without necessarily restrincting to the case of 
uniform distribution. The Hausdorff distance between two sets Ai and A 2 
is defined as 


jgn ^ 

\K\ 


dn(Ai,A 2 ) = inf{e > 0/Hi cd 2 + eBd{ 0,1) and A 2 £ A\ + eBd{ 0,1)}, 

where Bd(a, r ) is the closed d-dimensional Euclidean ball with center a e 
and radius r > 1. 

The Hausdorff distance between the random polytope K n and the sup¬ 
port K has been studied under various assumptions in the litterature. Barany 
[2] showed that if X\ ,..., X n are uniformly distributed in a d-dimensional 
convex body K with twice continuously differentiable boundary and positive 

Gaussian curvature, then E [d#(iF n , K )] is exactly of the order 

with constant factors that depend on K. In the planar case (d = 2), Ko¬ 
rostelev and Tsybakov [14] proved a similar result, with constants that do 
not depend on K , under the extra assumption that a ball of given radius 
rolls freely on the inner side of dK. This assumption is frequent in the 
litterature, see for instance [1,17] and [21] for more details about the free 
rolling ball property. Korostelev and Tsybakov [14] also proved that under 
the same assumptions, no random set achieves a strictly better uniform rate 


^lnn j 2/(d+1) 


2 



with respect to the expected 


of convergence than that of K n , i.e., ( — 

V n 

Hausdorff distance of K n to the true support K. 

Diimbgen and Walther [7] considered general distributions on convex 
bodies in order to show as particular cases that if the boundary of a convex 
body K satisfies some smoothness condition, then almost surely, 



d H (K n ,K) = 0^ l ™y +1 j 

if the random points are uniformly distributed in K and 

d H (K n ,K) = 0^^y' l) j 


if they are uniformly distributed on the boundary of K . In the case of 
uniform distribution in a convex body K , they also showed that under no 
smoothness assumption on the boundary of K , 




almost surely, 


with constant factors that depend on K and d. 

More precisely, when X \,..., X n are i.i.d. random points of OK with a 

smooth density and if dK is smooth enough, Glasauer and Schneider [10] 
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computed the limit in probability of 


-) dn^Kn, K). Schneider [19] 

n) 


conjectured that this limit holds almost surely. 

We propose uniform deviation inequalities for the Hausdorff distance 
between the random polytope K n and an enclosing convex body K of the 
support of the points Xi,..., X n under general assumptions on the behavior 
of that distribution near the boundary of K. In most cases, K will be equal 
to the support itself, or to its convex hull. By uniform, we mean that the 
bounds depend neither on K nor on the distribution itself. Particular cases 
include uniform probability measures in smooth convex bodies, in polytopes 
and on the boundary of smooth convex bodies. We also extend the opti¬ 
mality result proved in [14] in the case of uniform distributions to higher 
dimensions. 

This paper is organised as follows. After displaying our notation and 
preliminary definitions, we prove our main theorems in Section 2. The first 
one bounds the performance of the random polytope, under assumptions 
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made in the same spirit as in [7] and the second one concerns its functionals, 
including its mean and maximal widths. That section ends with a discussion 
on misspecification, i.e., about the non convex case. In Section 3, we derive 
from our main theorems upper bounds on the performance of the random 
polytope in three cases: uniform distribution in smooth convex bodies, in 
polytopes and on the boundary of smooth convex bodies. For the case of 
polytopes, we give non uniform results, in order to avoid making too complex 
assumptions on the facial structure. For that of uniform distribution in 
smooth convex bodies, we extend the optimality result of [14]. Section ?? 
is devoted to the proofs. 

1.2 Notation and introductory considerations 

Recall that throughout all this work, the dimension d is a fixed positive 
number. We denote by K.d the class of all convex bodies in R d and by k}J ^ 
the class of all convex bodies included in the closed unit Euclidean ball in 
R d . 

If p is a positive integer, the p-dimensional closed Euclidean ball with 
center a e MP and radius r > 0 is denoted by B p (a, r ) and the p- 1-dimensional 
Euclidan unit sphere is denoted by S' p_1 (so, 5 P_1 £ M p ). Th p-dimensional 
volume of B p { 0,1) is denoted by f3 p . 

We denote by || • ||, without any subscript, the Euclidean norm in R d . 

The interior of a set A is denoted by int(A). Its boundary is denoted 
by dA. 

The support function h,A of a compact set A c M. d is the mapping 
h'A • 

u nra x{(u,x)/x e A}, 

where (•,■) stands for the canonical inner product in R d . Note that if A is a 
compact set, then h a = h conv (A), where conv(A) is the convex hull of A. 

We recall the two following properties of support functions. If K is a 
convex body, then hx is positively homogenous, i.e., 

Vu e R d , VA > 0, Iik(\u) = XIik(u) (1) 

and it satisfies the triangle inequality: 

Vu, v e R d , Hk(u + v) < Hk(u) + Hk{v). (2) 

In particular, (2) implies the reverse triangle inequality 

Vu, v e M. d , hx{u - v) > Hk{u) - hx(v). (3) 
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We define the pseudo-norm || • \\k associated to a convex body K (also 
called the Minkowski functional of K ) as 

||a ;|k = nrin{A > 0/x e Ail'}, Vx e R d . 

This is a norm if and only if K is symmetric, i.e., K - -K. If I\ = S^ -1 , 
then || • || k is the Euclidean norm in 

The polar body K° of a convex body K is the convex set defined as 

K° = {x € M d /Vy e K, < x, y >< 1}. 

If 0 e int(K), then K° is compact, hence a convex body and 

Vx e K d , Hk(x) = ||x||ifo. (4) 

It is also true that for bounded K. 0 is an interior point of K°. 

Let us state a property of the Hausdorff distance, which will be of interest 
for our purposes. The Hausdorff distance between two convex bodies I\ and 
K' can be expressed in terms of their support functions: 

dii(K,K')= sup \h K (u)-h K '{u)\. 

If I\ e JCd, u e S d l and e > 0, we denote by Ck(u, s) the cap of K in the 
direction of u and of width e: 

Ck(u,£)-{x€K / (u,x) > hx(u)} . 

When there is no ambiguity, we may omit the subscript K , especially in the 
proofs. 

We denote by A(M rf ) the collection of all probability measures on 
equipped with its Borel cr-algebra. For // e A(M rf ), we denote by the 
corresponding expectation operator. The support of a probability measure 
fj, in is defined as follows. Let be the collection of all open subsets 
O £ with fi(O) = 0 and define = |^J O (set = 0 if O M is empty). 

Because is a metric space, this is the largest open subset of M. d mapped 
to zero by fi, inclusionwise. The support of fi, which we denote by Supp(fi), 
is the closed set defined as the complement of il lt . 

If a, L and £q are positive numbers with £q < 1, let 

M(a, L, £o) = {(/r, K) eA(M d ) x / Supp(p) c K, 

Vu e S^ 1 , Ve e [0, £q], p (Ck(u, e)) > Le a }. 
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Note that if (p,K) e A4(a,L,so) for some positive numbers a, L,£o> then 
K is necessarily the convex hull of supp(p). Indeed, K is convex and it is 
easy to see that every closed halfspace that contains either supp(p) of I\ 
necessarily contains the other as well. For the sake of simplicity, we may 
also refer to K as the support of p. 

Roughly speaking, if p is the uniform measure in a convex body K with 
smooth boundary, then the couple (p,K) belongs to the class M(a,L,e o) 

for a = ——— and L, £q suitably chosen. If p is uniform in a polytope K , then 
one can take a = d and if p is the uniform distribution on the boundary of a 
smooth convex body, a = ———. More details on this account will be given 
in the sequel. 

If n is a positive integer and X\,...,X n are i.i.d. random points uni¬ 
formly distributed in some convex body K , we denote by their joint law 
(for the sake of notation, we again prefer not to show its dependency on 
n) and by Ea; the corresponding expectation operator. Similarly, if the n 
i.i.d. random points are uniformly distributed on the boundary of a convex 
body K, we denote by Pqa' their joint law and by E gx the corresponding 


expectation operator. 

If (u n ) and (v n ) are two positive sequences, we write that u n = 0(v n ) 
or, equivalently, v n = u(u n ), if the sequence (u n /v n ) is bounded and u n x v n 


if u n = 0(v n ) and u n = co(v n ) simulatenously. 

Throughout the rest of the paper, we denote by C a = inf 

t> o 

any positive number a. 


(1 + t) a 
1 + i° 


> 0 for 


2 Estimation of the support of a measure 

2.1 Performance of the random polytope 

In this section, we quantify the performance of the random polytope K n 
in terms of deviation and moment inequalities for dn(K n , K). In other 
words, we ask how small this distance can be with high probability, or in 
expectation. We assume that the probability measure p of the i.i.d. observed 
points and the convex body K satisfy ( p,K ) e 7W(a, L, eo)> for some positive 
numbers a,L and £o- From a statistical point of view, the random polytope 
K n can be interpreted as an estimator of K, if K is unknown. In that 
case, one wish to recover K from the observed points X \,..., X n only. One 
advantage of the estimator K n is that it is adaptive to p, i.e., its computation 
does not require any knowledge about p. In Section 3.1.2, we will show that 
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under some precise condition, K n is an optimal estimator of K in a minimax 
sense. 


Theorem 1 . Let d and n be positive integers and a, L, <sq be positive numbers 


ith 0 < so < 1. Set n = max (1, —-V a n = ( — -) 

V C a aL) \ n ) 


with 


and b n = ri « . The 


random polytope K n satisfies the following uniform deviation inequality on 
the class A I(a,L, sq): For all nonnegative number x satisfying a n + b n x < £q 


sup n \dn{Kni K) ^ 2 a n + 2 b n x\ < 12 d exp ( -C a Lx a ) 


Remark 1 . • If a n + b n x > 1, the upper bound in Theorem 1 can be re¬ 

placed with 0, since for all ( n,K ) 6 M.(a, L,£q), 

dn{K n ,K)< 2 /x-almost surely. 


• For £q < a n + b n x < 1, the upper bound in Theorem 1 can be replaced 
with its value at x = (eo - Qnfbf 1 , since the left side of the inequality is 
a nonincreasing function of x. 

• It is important to notice that the constant factors in Theorem 1 are 
much smaller than in Theorem 1 in [6] when the dimension d becomes 
large. Unlike for those constants in [6], we believe that the factor 12 d 
cannot be much improved in general. 

Theorem 1 yields the following moment inequalities. 


Corollary 1 . Let d be a positive integer and a,L,£o be positive numbers 
with 0 < £q < 1. Then, for all real number q> l, 


sup [d H (K n , K) q ] =Ol ( —) 

U,K)aM{a,L,e 0 ) \ \ n / 


2.2 Estimation of functionals of K 

Previously, we focused on the performance of K n , as an estimator of K. The 
set K itself may not be of main interest to the statistician. In this section, 
we study a class of functionals of K and we bound the performances of the 
corresponding functionals of K n as plug-in estimators of those of K. In 
this context, a functional is a mapping T : /Q —> M. If K is unknown, 
a functional T(K ) can be estimated using the so called plug-in estimator 
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T{K n ). Plug-in estimators are often suboptimal. For instance, if T is the 
volume, i.e., T(K) = \K\ and p is the uniform probability measure in J\, [9] 
proved that the plug-in estimator \K n \ is suboptimal and proposed another 
estimator, based on a three-fold sample splitting, that achieves optimality 
in a minimax sense. Let us introduce some new notation. 

The unit sphere S^ 1 is equipped with the Borel er-algebra that it inherits 
from M. d . For any real number p > 1 and any measurable function / : S d ~ 1 -»■ 
M, denote by 

i 

\\f\\p = \f( u )\ Pda ( u )Y » 

where a is the surface area measure on the sphere with the normalisation 
a (S^ -1 ) = 1 and 

II /II oo = sup \f(u)\, 
provided these quantities are finite. 

If K e JCd, we still denote by hx the restriction of its support function 
to the unit sphere and we define 


<f>K ■ S d ~ l K 

u i-> hx(u) + hji(-u). 

We are interested in functionals of the form T p (K) = ||/i^|| p or S P (K ) = 
| cj )k ||p, for real numbers p > 1 or p = oo. Extensions to broader classes 
of functionals is possible. Note that S± is the mean width and S^ is the 
maximum width. Asymptotics of the mean width of K n are known in many 
cases. When p is uniform in a smooth convex body K, it is well known that 

— C K n^, 

n-* oo 

where ck is a positive number that depends on K , see [1,20] for explicit 
formulas. If p is uniform in a polytope K , Schneider [18] proved that 

5i(A')-E A '[5i(^ n )] — c'xn 1 *, 

J n^oo 

for another porisitve constant c' K that also depends on K. When p is sup¬ 
ported on the boundary of a smooth convex body K with a positive density 
g with respect to the surface area measure of dK , Muller [15] showed that 

Si(A')-E A -[Si(j? n )] —► c " K n 

L J n->- oo 



where the positive number c'g K depends on g and K. Here, we prove moment 
inequalities for S\(K n ) as well as for the other functionals S p ,T p ,p e [1, oo) u 
{oo}. In [1] it is shown that for the uniform distribution in a smooth convex 
body K, the variance of S\(K n ) is of the order n , up to constant factors 
that depend on K. Our results do not allow to recover such estimates of the 
variance because we do not bound moments of the centered versions of the 
functionals. 

Let pe[l,oo)u{oo} and q > 1 be a real number. If K and K' are two 
convex bodies, then by the triangle inequality, 


o- a, 

Ss 

1 

-e 

VI 

Cr 

E? 

i 

(5) 

\S P (K) - S p (K')\ q < \(f> K - (j) K ’\ q p . 

(6) 


In the next Theorem, we give uniform upper bounds for E^ [\\Jik - 
and E A( [\\(f>K ~ 4 >k lip] when (p,K) ranges in M.(a,L,e o) for given positive 
numbers a,L,e o. Note that the case p = oo is treated in Corollary 1 and we 
treat here only finite values of p. 


Theorem 2. Let d be a positive integer, a, L > 0, 0 < £o ^ 1 and p,q> 1 be 
real numbers. Then, 


and 


sup E„[|| h K -h R J g p ] 

(t-i,K)cM(a,L,£o) 


sup 

(p,A')eA4(a,L,£ 0 ) 


J| <t>K $K n 


Pi 

HpJ 


= o(n “ j 

= 0(n"«). 


As opposite to the case p - oo (see Corollary 1), the logarithmic factor 
can actually be dropped for finite values of p. We explain this phenomenon 
by the fact that L°° norms are more restrictive than L p norms (p e [1, oo)), 
since they capture the largest values of a function, no matter the size of the 
subset where large values are achieved. When p is the uniform probability 
measure in a smooth convex body, Theorem 5 in Section 3.1.2 shows that 
the logarithmic factor cannot be avoided for p = oo. 

By (5) and (6), next Corollary is a direct consequence of Theorem 2. 


Corollary 2. Let d be a positive integer, a,L > 0, 0 < £o ^ 1 and p, q > 1 be 
real numbers. Then, 

sup E /t [| T p {K) - T p (K n )\ q ] = O ini) 

(At,A')eA4(a,L,£ 0 ) 
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and 

sup E /t [|S p (A') - S p {K n ) |«] = O (ni\. 

( t i,K)eM(a,L,£ 0 ) V ' 

2.3 Misspecification: The non convex case 

In the setup of random polytopes, convexity of K is a natural assumption. 
However, in practice, this assumption is very strong. In particular, it could 
be misleading and the model misspecified: For instance, the probability 
measure of the observed points X \,..., X n may be supported in a non convex 
set or on the boundary of a non convex set. In that case, the random convex 
hull K n does not estimate the support K itself, but instead, its convex hull. 
Define the new classes M'(a,L,e o) as 

M'(a, L, e 0 ) = {(//, K) eA(R d ) xC* 1 ' / Supp(p) c K , 

Va e S' d_1 , Ve e [0,eq],h(Ck(u,e)) > Le° }, 

where is the class of all compact subsets of R d with positive volume. 
Then, all results stated in the previous sections (theorems 1, 2 and their 
subsequent corollaries) remain valid if Xi(a, L,eq) is replaced with the new 
class M.'(a, L,eq) and du(K ni K) is replaced with du{K n ,conv{K)). The 
key argument for this statement is the simple fact that the support function 
of a compact set is equal to that of its convex hull. This is of particular in¬ 
terest if one wishes to estimate functionals of a compact, but not necessarily 
convex set, such as its mean width. For instance, adaptation of Theorem 2 

to the non convex case shows that the mean width of a compact, non convex 

_ 

set can be recovered at the speed n “ under some assumptions on /a. A case 
of particular interest is when n is the uniform probability measure in a com¬ 
pact set K which has a freely rolling ball of a given radius r along the the 
inside part of its boundary (see Figure 1). In that case, similarly to what 

we will see in Section 3.1.1, the convex hull of K can be recovered by K n 

2 

( In n \ d + 1 

-1 and its functionals S P (K ) and T p (K), for p e [1, oo), 

_ 2 _ 

are estimated by the plug-in estimators at the speed n d +!. 
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Figure 1: The non convex case: K is a compact set with a free rolling ball 
on the inner side of its boundary, of give radius r. 


3 The case of uniform probability measures 


3.1 In a smooth convex body 

3.1.1 Performance of the random convex hull 


Let r be a positive number. We say that a compact set K c satisfies 
the r -rolling condition if at each point x of its boundary dK , there exists 
a Euclidean ball B with radius r, included in K and containing the point 
x (see Figure 1). For all real numbers r e (0,1], we denote by /C^ the 
class of all convex bodies that satisfy the r-rolling condition and that are 
included in the centered unit Euclidean ball. If K e /C^ and /i is the uniform 
probability measure in K, it is easy to see that (/ jl,K ) e A4(a, L,eq) for 


a = 


d+ 1 _ 2{3 d - 1 r d 2 1 


,L = 


Pd(d+ 1) 


and £q = r. Indeed, if B is any Euclidean ball with 
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radius r, then for all u e S d 1 and 0 < e < r, 


h(C k (u,£)) = 


I C K (u,e)\ 

\K\ 

I C B (u,e)\ 


f (x(2r-x)) 2 p^dx 

Jo 


1 re 

Jd 

Pd~ir C - 
Pd 

d -1 

2p d -ir 2 d±i 
-e 2 . 


X 


e d-1 

x 2 dx 


Pd(d+ 1 ) 

Therefore, the following result is a direct consequence of Theorem 1. 

d -1 

2p d -ir 2 


(7) 


Theorem 3. Let d, n be positive integers and 0 < r < 1. Let L = 


/ 


Pd(d+ 1) ? 


r = max 


1, 


dPd 


' t In n \ d+i _ 2 _ 

and b n = n d +\. The random 


/ r m n \ < 
\ n ) 


polytope K n satisfies the following uniform deviation inequality on the class 
KpJl: For all nonnegative number x satisfying a n + b n x < r 

sup F k \dfj(K n , K) >2a n + 2b n x\ < 12 fl! exp (^-Cd+i Lx a ^j . 


KtK 


As a consequence of this Theorem, we have the following uniform mo¬ 
ment inequalities. 


Corollary 3. Let d be a positive integer 0 < r < 1. Then, for all real number 

q> i, 

sup E K [d H (K n ,K) q ] =0 

a,r 

In the next section, we prove that these moment inequalities are tight in 
a minimax sense. 

As far as functionals are concerned, Theorem 2 can also be applied in 
the present case, yielding 
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Theorem 4. Let d be a positive integer and 0 < r < 1. Then, for all real 
number q > 1, 

sup E k [|| h K - hft n ||p] =o(n~d^^ 
kzkST* 

a,r 

and 

sup E k [\\(f> K - cj) Rn |||] = O (n~d+i ) . 

KzK. ( , 1) 

a,r 

This theorem shows that the plug-in estimators of the functionals T p (K) 
and S P (K), for real numbers p > 1, perform in expectation at a speed at least 
n~ 2 /( d+1 ), when K is a convex body with a free rolling ball of give radius. As 
we discussed in the beginning of the section, this result was already known 
for the mean width S±(K). 

3.1.2 Optimality of the random convex hull in smooth convex 
bodies 

In [5], we showed that the random polytope I\ n is rate-optimal in a minimax 
sense, with respect to the Nikodym distance. Let us recall briefly this result. 
If n is a positive integer and X \,..., X n are random variables mapping the 
probability space P) onto a given measurable space X , we call a set 

valued statistic any set of the form K n {X \,... ,X n ), where K n maps X n to 
the set of all compact subsets of M rf , such that the function 

HxM 11 -. {0,1} 

1 ifx€K n (X 1 (u),...,X n (u)) 

0 otherwise 

is measurable with respect to the product of T and the Borel u-algebra of 
M. d . For the sake of notation, we write K n instead of K n (X\,... ,X n ). For 
instance, if X = and X\,... , X n are i.i.d. random points, then K n is a 
set valued statistic. If a set valued statistic K n is aimed to approximate K, 
which may be unknown in practice, it is called an estimator of K and the 
variables \K n A K\ and dn{K n ,K) are (measurable) random variables. A 
more general definition of set valued random variables can be found in [16], 
but it is not needed for our purposes. In [5], it is shown that for any 

never converges to zero faster than 

when K n = K n , when n -> oo. The random polytope K n is called rate- 
minimax optimal on the class /Q, for the rescaled Nikodym distance. Next 


estimator K n of K, sup E# 
KeK. d 


| K n A K | 
1*1 
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theorem, together with Corollary 1, establishes rate minimax optimality of 
K n on the class /c[, 'J, for the Hausdorff distance djj- 

Theorem 5. Let d be a positive integer and 0 < r < 1. Then, for all real 
number q > 1 and any set estimator K n , 

sup E k [, d H (K n , K) q ] = u 

K ^ {1 1 

r,a 

We restrict r to be strictly less than 1 in Theorem 5, since when r = 1, 
the class /Cj 1 j is made of only one element, B^( 0,1) and the trivial estimator 
K n = Bd( 0,1) satisfies sup E# [dn(K n , K )] = 0. 

In [5] and [6], the risk of K n was measured using the Nikodym distance 

I K n A K\ 

rescaled by the volume of K , namely, ——-. Rescaling with the volume 

of K allowed to consider convex bodies K not necessarily included in some 
given bounded set such as ^(0,1). Same is possible in our present case 
if the definition of the Hausdorff distance is slightly modified. This is the 
purpose of next section. 



3.1.3 Rescaled performance of the random convex hull 

For the uniform case, we propose another way of measuring the distance 
between the support K and K n , which allows to remove the restriction 
K c Bd(0, 1). Namely, we define the distance between two subsets K 
and K' of as 

d L (K,K') = 

inf{e > 0/3a e R d , a + (1 - e)(K - a) c K' c a + (1 + e)(K - a)}. 

Note that does not properly define a distance since it is not symmetric. 
However, if di,{K,K') is sufficiently small, then dL(K,K r ) and dr{K',K) 
are roughly the same. One advantage of dL over dn is that it is invariant 
under any invertible affine transormation T : K') = di(T(/i), T(K')). 

In turn, considering dL instead of du allows to get rid of the assumption 
K c Bd{ 0,1) in our results. 

For 0 < r < 1, let us denote by K- r ^ the class of all convex bodies in 

\ 

-] . By introducing the 

Pd J 
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number T in this definition, we ensure that all convex bodies that are similar 
to K , i.e., obtained by rescaling, translating, rotating or reflecting K, will 
be in /Q. r as soon as K is. 

Theorem 6. Let d and n be positive integers and r e (0,1]. Set the posi¬ 


tive constant 73 = max 1 , —-—- 1 — , a n = 

\ p d -iCd±ir d L J 

Then, the random polytope K n satisfies 


sup P k \dL(K n , K) > 2 a n + 2 b n x\ < 6 d exp (- 73 .T 2 
KeK r , d v 


(*r) 


2 

, d +1 


and b n = n 


2 

' d+l 


for all nonnegative real number x with a n + b n x < 1. 

Remark 2. For K e KL r ^, since K n £ K Fk- almost surely, di,{K n ,K ) < 1. 
Hence, the right side of the inequality in Theorem 6 can be replaced with 0 
if a n + b n x exceeds 1. 

Theorem 6 , together with the previous remark, implies the following 
moment inequalities. 

Corollary 4. Let d be a positive integer and r e (0,1]. Then, for all real 
number q> 1, 

sup E K [d L (K n ,Ky]=o(( — Y +1 ). 

KeK r , d \\ n ) J 


3.2 In a polytope 

The polytopal case is more intricate and we choose only to present the 
heuristics here. Let K be a polytope and u e S d ~ l . Let H be the (unique) 
supporting hyperplane of K in the direction u. Then, K n H is a face of 
K and let k be its dimension, i.e., the dimension of its affine hull. As 
e -» 0, Pa' [Ca'(u,£)] goes to zero as a polynomial of degree d- k in e. This 
speed of convergence is the fastest when k = 0, i.e., when H supports K 
at a vertex. By a compacity argument, it follows that for all polytope K , 
(P k,K) e Ai(a, L,£q) with a = d and L and £0 being two positive numbers 
depending on K. Hence, 

E K [d H (K n ,K)] = 0^^ 1/d y 
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with constant factors that depend on K through the parameters L and So- It 
is possible to write uniform deviation and moment inequalities by restricting 
the class of polytopes in such a way that L and £q do not depend on K 
anymore, but the corresponding assumptions (similar to the free-rolling ball 
condition in the cas of smooth convex bodies) would be quite technical. 

One comment seems necessary. In the case of uniform observations, 


Groemer [11] showed that when K ranges in K. d , Ex 


j Kn A K | 
1*1 


is maximal 


when K is an ellipsoid (any ellipsoid, by affine invariance) and it was con¬ 
jectured by Barany and Buchta [3] that it is minimal when K is a simplex. 
What is known (see [4]) is that it converges to zero at the fastest rate (in 
terms of n) when K is a polytope. However, our results suggest that the 
opposite is true for EA'[d^(* n , K)]: we showed that it converges at a speed 
Tnn\ 2/(d+1) , „. , /In r0 1/d 


at least 




when K is a Euclidean ball and at at least 




when K is a polytope. Hence, for the Hausdorff distance, the situation seems 
to be reversed as compared with the Nikodym distance. 


3.3 On the boundary of a smooth convex body 

In this section, we assume that d > 2 and we consider convex bodies in 
that satisfy a rolling condition. Let us fix a positive number r e (0,1] and 
let K e /C^*J. We consider n i.i.d. random points X\, ..., X n uniformly 
distributed on dK, where n is a positive integer. Let u e S^ -1 and 0 < e < r. 
Then, the probability content of Cx(u,e), with respect to Pax, is the ratio 
between the surface area of dK n Ck(u, e) and that of dK. By the r-rolling 
condition, the surface area of dK n Ck(u ,£) is greater or equal to that of 
a spherical cap of radius r and height e. If we denote by Ad the surface 
area of the unit sphere the surface area A(r, e) of such a spherical cap 
satisfies 

A(r,e) = \A d r d - 1 £^ t^(l- ty^dt 

>-A d r d - 1 n&dt 

2 Jo 

d -1 d -1 

> A d r 2 £ | 

In addition, since K is convex and is included in B d ( 0,1), the surface area 
of its boundary is bounded from above by that of S d ~ 1 , which we denoted 
by A d . In turn, 

P dK [C K (u,s)] > 
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and therefore, (P qk,K) belongs to the class M. ^—-—,r 2 ,rj. Hence, as a 

consequence of Theorem 1, we have the following result for the convex hull 
of independent random points uniformly distributed on the boundary of a 
smooth convex body: 


Theorem 7. Let d > 2, n be positive integers and 0 < r < 1. Define 

/ \ - 2 


t 4 = max 


1 , 


2d 


^ Ca -1 (d - l)r 2 j 


r 4 In n \ d-i ^ b - & _ The mn _ 




dom polytope K n satisfies the following uniform deviation inequality on the 
class : For all nonnegative number x satisfying a n + b n x < r 


sup P dK \dn(K n ,K) > 2a n + 2b n x] < 12 d exp r 2 x a j . 


KeIC 


(!) 

d,r 


Uniform moment inequalities follow from this theorem: 


Corollary 5. Let d > 2 be a positive integer and r e (0,1]. Then, for all 
real number q > 1 , 


sup E dK [d H (K n , K) q ] 
( 1 ) 


O 


K^K 


r,d 


((^) 


2q_ 


Remark 3. Section 3 is devoted to uniform probability measures, either in 
a convex body K, or on its boundary. In both cases, the uniform probability 
measure admits a constant density f, with respect to either the Lebesgue 
measure or the surface area measure of dK. From the proofs, it is clear that 
our results extend to the case when 0 < a < f(x) < b for all x in the support 
of /i (either K, or dK), where a and b are two fixed positive numbers. The 
rates of convergence, in terms ofn, are still the same and only the constants 
are modified, depending on a and b. 


4 Proofs 

We start with the following lemma, which is directly inspired by Lemma 5.2 
in [ 8 ]. 

Lemma 1. Let K e /Q and 5 be a positive real number. There exists a 
subset AIs of dK, satisfying the following: 
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/ 3 \ d 

• The cardinality of Ms is not larger than I -1 ; 

• For all u e 8K, there exist two sequences (uj)j >o £ Ms and (Sj)j> 1 £ M 
such that the following holds: 

— Vj > 1,0 < Sj < 5i and 

oo 

— u = u 0 -Y / djUj. 

3 = 1 

Proof of Theorem 1 Let (//, Lf) e Ml(a, L, £o). Let z € [0,£o] and 6 = e/4. 
For all u € the following statements are equivalent: 

• hg n (u) < h K (u) - e 

• All the points X±,.... X n are outside Ck(u,z). 

Hence, 


n[h Rn (u) < h K {u)-e] = At [Vi = 1,..., n, X. t { C K (u,z )] 

= (l- M (Cx(u,e))) n 

<exp(-Lne a ). (8) 

Let Ms be a subset of S rf_1 satisfying Lemma 1, for the unit Euclidean ball. 
Let us denote by A the event {Vu e Ms,h^ n (u) > hx^u) - e}. By (8) and 
the union bound, 

/x(A)>l-Qj exp(-Lne a ). (9) 

OO 

Let A hold and u e S^ -1 . By Lemma 1, we can write w = tto + ^ S-j Uj where 

j=i 

( u j)j> o - -Ms and 0 < Sj < 8 3 , Vj > 1. Note that /j- almost surely, K n 9 K c 
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Bd( 0,1), hence hj> n (u) < hx(u) < 1 ,Vu € S d 1 . Thus, 


^ l K n ( u ) _ h Kn ( u o + E &j u j) 

j =i 

oo 

^ h R n ( u o) - E 5 ]}i kS~ u ^ b y (i) and ( 2 ) 

i=i 

> h K (u 0 ) -£- r 

1-5 

°° • 5 

> 1ik(u) - E fi J hK(uj) ~ £ -- again by (1) and (2) 

i=i 1 “ ^ 

, , x 25 

> riK\u) - e ■ 


5-1 


> hx(u) - 2s, 

since 5 = e/4 < eo/4 < 1/2. 
Hence, by (9) and (10), 


( 10 ) 


/J [djf(if n , ^0 > 2e] = n [3u e 5 d 1 ,hg n (u) < h K (u) - e] 

< (y) d exp(-Lne“) 

< 12 rf exp (-Lne Q - dine). 


( 11 ) 


By setting e = a n + b n x for some nonnegative number x, where a n and b n are 
defined in Theorem 1, we have that 


£ a >C a 


( T\ In n x a \ 

-+ — h 

V n n ) 


and 


In e > In a n > - 


Inn 


a 


The conclusion of Theorem 1 follows. 


□ 


Proof of Corollary 1 Let Z be a nonnegative random variable and q be 
a positive number. Then, by Fubini’s theorem, 

Jf' OO 

t q ^F[Z > t]dt. 
o 
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Let us apply this equality to Z = du^Kn-, K) (which we denote by dn for 
the sake of simplicity) and q> 1. Since dff(K n , K) < 2 /i-almost surely, 

E »[d H {K n ,KY] 

X oo 

[dn > t] dt 

J f' 2 dn r 2eo 

t q ~ fi [dn >t]dt + q t. q ~ /x [dn > t] dt 

0 J2a n 

+ q f tJ~ l n[dn >t]dt. (12) 

J2£n 


In the first integral, let us bound the probability by 1. In the second and 
third integrals, we perform the change of variable t = 2 a n + 2 b n x. Then, by 
Theorem 1 and the remark that follows it, (17) becomes: 


E^[d H (K n ,K) q ] 


rCo-anlL, 1 

(2a n ) q + 2qb n J (2a n + 2b n x) q n[du >2a n + 2b n x]dx 


2qb n f 


(1 -a n )b„ 


(eo-a^b^ 

(£o~a n )b~ 1 


(2 a n + 2b n x) q 1 /i [dn > ^dn + 2 b n x] dx 


r (£ 0 —a n )O n 

< (2 a n ) q + qb n / (2a n + 2b n x) q x 12 d exp (-C a Lx a ) dx 

Jo 

+ q2 q ~ 1 11 - e 0 |6)) 1 12 d exp (~C a L\l - e 0 \ a b~ a ). 

(13) 

In the second term, we bound (a n + 6 n x) 9_1 by 2 q ~ 1 a^ 1 + 6^ _1 x 9-1 , which 
yields 


X 


(eo -a n )b n 


(2 a n + 2b n x) q 1 12 rf exp (- C a Lx a ) dx 
< 12 d A q ~ 1 exp ( -C a Lx a ) dx + 6^ 1 x q ~ l exp ( -C a Lx a ) dxj . 


(14) 


Since b n = 0(a n ), altogether, (18) and (14) yield Corollary 5. 


□ 


Proof of Theorem 2 Let e M(a,L,£o) and let u e S d 1 . By (8), 

for all e e [0 ,£q], 


H \hft (u) < h K (u) - e] < exp (-Lns a ). 
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For e e (e 0 ,2], 

/j < hx(u) - e] < exp (-LtiEq ), 

since the left hand side is a nonincreasing function of £ and if e > 2, 

M [%„(“) ^ M«) -e] = 0. 

Hence, by a similar argument as in the proof of Corollary 1, for all real 
number k > 1, there is some positive number that depends neither on 
/i, K nor on u and such that 

E/u [|h*r(u) - ^ n ( u )! fe ] - c fc n_ “ • (15) 

£ 

If g < p, Jansen’s inequality applied to the concave function [0, oo) 
and the expectation operator E^ yields 

[II h R n ~ h x\\p\ = [f sd ^ \ h k(u) - h Rn (u)\ p da(u ) j 
r r ~\ q / p 

^^[J sd _i\ hK (u)- h Rn (u)\ Pd <T(u)\ . (16) 

By Fubini’s theorem, which allows switching the expectation and the integral 
and by (15) with k = p, (16) yields 

which proves the theorem in the case q < p. If now q > p, the mapping 

£ 

[0, oo) b x h* xp is convex and Jensen’s inequality applied to the operator 
f S d- 1 yields 

(f sd J hK ^- h Rn^\ P ) * f sd J h K^)-h Rn ^)\ q , 

so again by Fubini’s theorem and using (15) with k = q, 

[ll^j? n “ lip] ^ c q n “ > 

which ends the proof of the theorem. □ 
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Proof of Theorem 5 We generalise the method used in [14] in the planar 
case to any dimension. This is a standard method that was already exploited 
in [13]. We fix R e (r, 1] and set Gq = Bd{ 0, R). Let 6 be a positive and small 
enough real number. For u e S d ~ 1 , we define the set G(u ) as follows. Let 
r ]: M -> M be a given nonnegative function, twice continuously differentiable, 
such that 


• Vx € M, |x| > 1 => r/(x) = 0 and 

• maxr/(x) = rj( 0) = 1. 

ireR 

For instance, one can take ?y(x) = e 4 g(2x - l)g(2 - 2x). Vx e M, where 


g:R-> 


X 


exp(-l/x) if x > 0 
0 otherwise. 


Let H be the supporting hyperplane of G at the point Ru e 8Gq ■ Identify this 
hyperplane with with origin at the point Ru. Then, a parametrisation 
of Go is 

G 0 = {(t,y)€HxR / \t\<R, R- \J R 2 - |f| 2 <y<R+\fR 2 - |t| 2 } , 

where we denote by j ■ j the Euclidean norm in H. We define the set G(u) 
by modifying the parametrisation of <9 Gq in a small neighborhood of Ru: 


G{u) = {(t,y) € H xf / \t\<R, 

R-^R 2 - \t\ 2 + ad 2 r, (^) <y<R + ^R 2 - |t| 2 }, 

where a is a positive number that can be tuned inependently of 5 so that 
G(u) e JCrrf, for small enough 5. Note that for all u e S d ~ l ,G(u) 9 Go and 


\G 0 \G(u)\ = a5■ 


r 41 )^ 

JR d \Rd) 


a6 d + l R d~l 


= cS 


2 d ~i 
d +1 


L ,mdt 


(17) 


for some positive constant c that depends on R and d only. In addition, if 
( u,u') e S d ~ 1 x 5 d_1 is such that ||w - v !|| > S, 

d ff (G(u),G(u')) = aS 2 . (18) 
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Let now J\f be a maximal 5-packing of S d 1 . i.e., a finite subset of S d 1 with 
maximal cardinality, satisfying 

V(it,u) e JV, ||it - u|| < 6 => u = v. 

Let us denote by N its cardinality and by u \,... ,un its elements and, for 
j = 1,.. ., N, set Gj = G(uj). By a standard argument (e.g. [12]), N > c 5 ~ d , 

i 

( Inn\ d+i 

-1 , we apply Lemma 

n / 

1 in [13] to the sets Go, Gi,..., Gat, whose assumptions are satisfied thanks 
to (17) and (18) and this proves Theorem 5. □ 


Proof of Theorem 6 The proof of this theorem follows the same lines 
at that of Theorem 1. Let d > 1, r e (0,1] and K e !C r ( j. For notational 
convenience, we set 



Let n be a positive integer and X, X i, X 2 ,..., X n be n +1 i.i.d. random points 
uniformly distributed in K. Without loss of generality, one can assume that 
0 is an interior point of K and f?d(0,r) £ K. Indeed, since K e there 
exists a point a e K such that f?d(a,r) £ K. One can then replace K 
with K - a and not change the distribution of du{K n ,K). Hence, for all 
u € S d ^ 1 ,hx{u) > r and 

d L (K n , K) < inf{e > 0/(1 - e)K c K n c (1 + e)K} 

= inf{e > 0/(1 - e)K c K n } 

= inf{c > 0/Vtt e S' 1-1 , (1 - e)liK(u) £ hj> n (u)}. 

Let e € (0,1] and u e S d ~ 1 . By the F-rolling condition, we have, similarly 
to (7), 

r er , - d -1 

|Gk(u,£?)I = \/x(2r - x) fid -1 dx 

a -1 r £r 1 

> /3rf_ir 2 / x 2 dx 

Jo 

2/3 d - 1 r d d ±1 

> —--e 2 

d + 1 

_2^_ 1 r d |iF| fei 

" /3d(d+l) ^ 2 ' 
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Hence, if r / = 


r d 

Md+iy 


h K,S U ) 

h K {u ) 


<l-£ 


= P [hg n (u) < h K (u) - eh,K{u )] 

< P [h Rn ( u ) < h K (u) - er] 

/ |C^,sr)| \ w 

" l 1*1 ) 


< exp 


1*1 

( d±l\ 

I -rjne 2 l. 


It follows, by positive homogeneity of support functions (see (1)), that for 
all u e M d \{0}, 


h K,S u ) 

h K (u) 


<l-£ 


< exp 


I d+l\ 

I -?yne 2 1. 


(19) 


We now introduce a finite subset of d(K °) satisfying the requirements 
of Lemma 1, then apply (19) to each of its elements and use a union bound 
inequality. By an analytic argument, we will then be able to bound from 

f h'R {u) ) 

above the probability that the inequalities I — n < 1 - e 1 hold simulta- 

[ h K (u) J 

neously for all u e d(K °) and, as a consequence of positive homogeneity of 
support functions, for all u e S d ~ l . 

Set 5 = e/2 and let Ms be a subset of d(K °) satisfying the requirements 
of Lemma 1. Applying (19) to each element of Mg and using a union bound 

/ 6 

inequality, yield that with probability at least 1 - I - 



hr> ( u) 

VutMe, n : ! > 1 


( 20 ) 


h K (u) 

Let us assume that Event (20) holds and let u be any element of d(K°). 

OO 

Write u= uq- £jUj, as in Lemma 1. Using (1), (2) and (3), 

3 =1 


^ l K n ( U ') ~ ^K n ( U 0 ) X! y^ l K r S U j)' (21) 

3 = 1 

In addition, for all nonnegative integer j, h,K(uj){l-e) < h^(uj) < hx{uj). 
Since 0 e int(K) and Uj e d(K°), hx(uj) = !%||k° = 1 and the following is 
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true, for all u e d(K°): 


h R n ( u ) > 1-e - 

3 = 1 

> 1 - e->l-2e. (22) 

1-5 77 

Finally, note that 0 e int(K° ), for K is bounded. So, for all u e S'^ 1 , there 
is a positive A for which Xu e d(K°) and 

hf> (u) hf> (Xu) 

1ik(u) hx(Xu) 


Hence, we showed that with probability at least 1 - 



) 


, , hf> (it) 

VueS^ 1 , , Kn ; . > l-2e. 

h K (u) 

It follows that 

dL(K n , K) < 2 £. 

Hence, we showed that 

P[d L (i? n ,K)>2e]<6 d exp (-Tyne 2 -dlnej, (23) 

which yields Theorem 6 by setting £ = a n + b n x. □ 
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